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粗糙 表面 接触 的 层级 模型 
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摘 要 :真实 的 材料 加 工 表 面 的 粗糙 度 往 往 跨 越 了 多 度 ar le ee dis 
级 接触 模型 。 对 于 较 大 尺度 部 分 ry 各 .公有 介质 力学 的 方法 计算 
相应 的 近似 接触 面积 和 局 部 压力 的 分 布 ;对 于 较 小 尺度 部 分 ,本 研究 基于 分 形 理论 ,利用 二 维 Can- 


tor 集 来 描述 层 晤 的 短波 长 微 西 峰 , 然 后 对 在 较 大 尺度 ss 他称 十 ， 计算 

?了 最 终 的 接触 面积 。 为 了 验证 所 建立 的 模型 ,将 得 到 的 整体 接触 响应 与 有 限 元 法 计算 结果 进行 了 
比较。 结果 表明 ,该 模型 能 够 充分 预测 粗糙 表面 的 弹性 接触 响应 ,而 且 大 大 降低 了 多 尺度 粗糙 表面 
接触 力学 建 模 的 计算 规模 。 这 对 进一步 研究 跨 尺度 粗糙 表面 之 间 的 接触 和 摩擦 具有 重要 意义 。 
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A hierarchical contact model for rough surface 


YAN Duo ,LI Leitao, LIANG Xuanming, XING Yuzhe 


(Department of Engineering Mechanics,SVL ,Xian Jiaotong University ,Xiran 710049 ,China) 


Y Abstract :Real machined surfaces are always rough across a wide range of length scales. This paper devel- 

ops a hierarchical contact model ,in which surface power spectrum density function is employed to divide 
roughness into two parts. In the macroscopic part, by cutting off short-wavelength roughness ,the approxi- 
mate contact area and the distribution of normal force are obtained adopting the continuum mechanics solu- 
tion. In the microscopic part ,the superposition of asperities for short wavelength is transformed into a two 
dimensional Cantor set based on fractal theory ,and the approximate contact area is further refined based on 
the local pressure and geometric parameters. The contact response of the whole surface is thus obtained by 
adding the modification of the microscopic part into the macroscopic part. For verification ,FEM is adopted 
to calculate the response of the fractal rough surfaces. The results demonstrate that our model can ade- 
quately predict the contact response of the multi-scale rough surface. The current study reduces the compu- 
ting scale of the mechanical model of rough surfaces and thus provides a new method modeling rough sur- 
faces in a larger scale span. This may have significant implications for further research on contact and fric- 
tion between multi-scale rough surfaces. 
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第 2 其 关 铎 ， 等 : 粗 烽 

大 量 的 实验 表明 , 当 2 个 材料 表面 发 生 接触 
时 ,处 于 真实 接触 状态 的 区 域 只 占 名 义 接触 面积 的 
一 小 部 分 。 此 外 ,在 微米 尺度 和 纳米 尺度 上 都 可 以 
观察 到 表面 粗糙 不 平 的 现象 。 随 着 观测 尺度 的 细 
化 ,许多 粗糙 表面 的 结构 表现 出 了 显著 的 分 形 特 
征 呆 。 预 测 跨 尺度 粗糙 表面 的 接触 响应 是 接触 和 摩 
擦 领域 的 一 项 具有 挑战 性 的 任务 。 

在 过 去 的 几 十 年 里 ,学 者 们 已 发 展 了 许多 描述 粗 
糖 表 面 形 貌 和 接触 响应 的 理论 模型 。ARCHARD'* 
将 粗糙 表 面 建 模 为 半球 释 加 在 较 小 尺度 上 的 自 相似 
结构 。GREENWOOD 等 “将 粗糙 表面 描述 为 曲率 
相同 且 高 度 服从 Gauss 统计 分 布 的 球面 。 此 外 ， 
NAYAKC 将 粗糙 表面 描述 为 二 维 随机 过 程 ,并 推导 
出 工 高 斯 型 随机 粗糙 表面 的 统计 特征 。 基 于 该 随机 
过 程 理论 ,BUSH 等 ”同时 考虑 了 粗糙 峰 高 度 和 曲 
举 揭 统计 特征 。 然 而 ,这 些 多 微 凸 模型 中 的 参数 是 
其 有 尺度 依赖 性 的 ,对 于 不 同 的 观测 尺度 不 具有 唯 
-他 ” 。 
一 基于 分 形 理论 ,MAJUMDAR 等 采用 W - M 函数 
来 播 述 粗糙 表面 ,并 得 到 了 整体 接触 响应 与 分 形 参 
效 击 间 的 关系 "1 。YAN 等 "进一步 建立 了 三 维 粗 
糙 尖 面 的 分 形 接触 模型 。 金 守 峰 等 "2 在 此 基础 上 
加 了 弹 塑 性 变形 的 修正 并 建立 了 修正 模型 。 赵 永 
武 区 "依据 微 凸 体 变形 的 连续 性 和 光滑 性 原理 提 
出 了 新 的 接触 模型 。 甘 立 等 ”通过 推导 第 一 、 第 二 
临 昭 弹 塑 性 接触 面积 与 临界 弹性 接触 面积 的 关系 ， 
建委 了 修正 的 弹 塑性 接触 模型 。 
- 亿 功 率 谱 常 被 用 来 分 析 粗 糙 表 面 的 几何 特征 。 自 
仿 射 分 形 粗糙 表面 接触 模型 中 ,采用 了 满足 寡 律 关 
系 的 功率 谱 " 。 然 而 ,实验 测量 的 结果 表明 粗糙 表 
面 的 功率 谱 函 数 与 尺度 的 竹 律 关系 在 较 大 尺度 下 并 
不 成 立 "" 。HYUN 等 ”发 现 ,引入 大 波长 截止 将 导 
致 真 实 接触 面积 与 载荷 的 比例 系数 的 小 幅 增加 。 因 
此 ,假设 粗粮 表面 在 所 有 长 度 尺度 下 的 形 貌 都 是 分 
形 的 ,这 可 能 是 不 合适 的 。 

数值 模拟 是 研究 粗糙 表面 接触 问题 的 有 力 手 
段 !"” 。 在 已 有 的 数值 计算 工作 中 ,考虑 的 表面 粗 
糙 度 仅 限 于 较 窗 的 尺度 跨度 。 在 有 限 的 计算 资源 
下 , 随 着 粗糙 度 的 细 化 ,可 以 模拟 的 模型 尺寸 将 急剧 
缩小 。 鉴 于 此 ,有 人 提出 了 解决 这 个 问题 的 方案 。 
例如 ,AN 等 中 开发 了 一 种 光谱 插值 方法 来 细 化 表 
面 。 然 而 ,这 种 插值 方法 无 法 在 小 尺度 上 提供 充足 
的 信息 。 多 尺度 方法 是 克服 这 一 困难 的 一 种 合理 有 


效 的 途径 。WANG 等 ” 建立 了 一 个 多 尺度 接触 模 
型 来 模拟 光滑 球体 与 实际 加 工 粗糙 表面 之 间 的 弹 塑 
性 接触 ,其 中 表面 总 位 移 场 由 全 局 水 平 的 体 变形 和 
子 级 的 球形 微 凸 体 变 形 组 成 。 

在 本 研究 中 ,基于 粗糙 表面 的 功率 谱 分 析 ,提出 
了 一 种 新 的 方法 来 解决 多 尺度 粗糙 表面 的 接触 问 
题 。 利 用 这 个 方法 ,计算 了 各 向 同性 粗糙 表面 与 刚 
平面 接触 时 的 法 向 载荷 与 真实 接触 面积 之 间 的 关 
系 。 结 果 表 明 ,用 此 方法 所 得 结果 与 直接 有 限 元 数 
值 计算 的 结果 吻合 较 好 , 且 大 大 提高 了 计算 效率 。 
这 将 有 助 于 进一步 研究 粗糙 表面 在 更 宽 的 频率 跨度 
内 的 接触 响应 。 


1 模型 描述 


实际 粗糙 表面 包含 多 尺度 粗糙 度 。 将 粗糙 度 看 
作 是 不 同 波长 正弦 波 的 到 加 ,由 粗糙 表面 的 高 度 信 
息 可 以 计算 得 到 其 功率 谱 密 度 ( power spectral densi- 
ty ,PSD ) 函数 。 对 于 典型 的 理想 各 向 同性 粗糙 表面 
的 PSD 函数 , 当 频 率 小 于 临界 频率 g, 时 ,PSD 函数 几 
乎 保持 一 个 恒定 的 值 P,。。 这 描述 了 由 于 加 工 轮 廓 
和 表面 处 理 的 影响 导致 粗糙 表面 在 较 大 尺度 上 呈现 
出 的 非 分 形 特征 。 当 频率 大 于 临界 值 9, 时 ,PSD 机 
数 将 与 频率 呈 对 数 线性 的 关系 ,这 意味 着 表面 在 细 
观 尺 度 上 呈现 分 形 特征 。 由 于 在 实际 粗糙 表面 的 测 
量 中 ,实验 采样 的 尺寸 是 有 限 的 ,而 采样 间隔 决定 了 
横向 的 分 辨 率 , 表 面 粗糙 的 频率 9 被 限制 在 qo 到 gy， 
的 范围 内 。 本 研究 所 提 到 的 频率 即 指 的 是 描述 空间 
频率 的 波 数 。 一 般 来 说 ,粗糙 表面 的 PSD 函数 可 以 
描述 为 


Po q. <qg 
2-2H 
C(q) = P, (4) q, <q<g (1) 
0 else 


其 中 ,五 为 粗糙 表面 的 Hurst 指数 ,可 以 通过 拟 合 被 
测 面 的 PSD 函数 得 到 。 

基于 上 述 PSD 函数 的 特性 ,本 研究 提出 一 个 层 
级 模型 来 计算 实际 接触 面积 ,如 图 1 所 示 。 以 某 一 
频率 q. 为 界 将 表面 粗糙 度 划分 为 大 尺度 和 小 尺度 两 
部 分 。 对 于 式 (1) 这样 的 PSD 函数 ,临界 频率 9, 是 
一 个 直观 的 选择 , 即 令 gq, = 9,。 首 先 ,使 用 仪 包含 大 
尺度 的 粗糙 峰 的 近似 粗粮 表面 ( 即 过 滤 掉 频率 高 于 
4. 的 粗糙 度 ) 来 估计 局 部 接触 压力 的 分 布 和 近似 接 


投稿 网 站 :http://ejam. xjtu. edu. en 微 信 公 众 号 :应 用 力学 学 报 


336 应 用 力学 学 


触 面 积 。 然 后 ,利用 基于 分 形 理 论 的 二 维 Cantor 集 
描述 小 尺度 部 分 的 粗糙 峰 , 并 对 其 在 接触 中 的 影响 
进行 了 分 析 , 对 大 尺度 部 分 得 到 近似 接触 面积 进行 
了 细 化 ， 0 全 出 了 考虑 多 尺度 粗糙 度 的 真实 接触 
面积 。 将 详细 描述 这 两 部 分 的 计算 过 程 。 


全 


频率 g/m 


Cm 


Cm 


gq, 


频率 q/m 


频率 g/m 


图 1 基于 粗糙 表面 PSD 函数 的 层级 模型 
Cig. 1 Hierarchical model based on the PSD function 
AN 


1QD 大 尺度 部 分 
ON 


二 在 这 部 分 ,包含 全 部 粗糙 度 信息 的 原始 粗糙 表 
将 通过 截 去 频率 高 于 9 的 粗糙 峰 的 方式 进行 重 
构 加 首先 ,通过 二 维 传 里 叶 变换 ,可 以 得 到 原始 粗糙 
表面 的 PSD 函数 和 临界 频率 9,。 然 后 ,将 频率 高 于 
gq. 的 各 频率 粗粮 的 振幅 置 零 ,而 频率 低 于 9 的 振幅 
与 贿 位 保持 不 变 。 最 后 ,对 截断 的 PSD 函数 执行 傅 
里 叶 逆 变换 得 到 仅 包 含 较 大 尺度 粗糙 峰 的 粗糙 表 
面 。 重 构 的 表面 只 包含 了 较 窗 的 频率 跨度 的 粗糙 ， 
并 且 细 观 的 形 貌 变 得 光滑 连续 。 因 此 ,可 以 使 用 有 
限 元 的 方法 对 重 构 表 面 进行 建 模 ,并 快速 计算 出 其 
与 刚性 平面 接触 响应 ,得 到 接触 过 程 中 的 法 向 压力 
分 布 和 该 有 限 尺度 下 的 近似 接触 面积 。 

进一步 假设 ,在 仅 包含 大 尺度 粗糙 的 接触 过 程 
中 得 到 的 局 部 法 向 压力 分 布 与 原始 粗糙 表面 接触 的 
保持 一 致 。 需 要 注意 的 是 ,本 小 节 得 到 局 部 法 向 压 
力 与 接触 区 域 的 最 小 尺度 是 在 有 限 元 网 格 的 尺度 ， 
离散 间距 就 是 粗糙 表面 网 格 的 采样 间距 , 即 V2n/g,。 
通过 有 限 元 计算 ,可 以 获得 每 个 处 于 接触 中 的 离散 
面积 单元 与 其 上 承受 的 法 向 压力 。 然 后 ,每 个 接触 
单元 上 的 真实 接触 面积 将 由 下 一 节 考虑 更 小 尺度 粗 
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粹 特征 的 计算 得 到 。 最 后 ,将 每 个 单元 上 的 真实 接 
触 面 积 相 加 , 即 可 得 到 在 一 定 法 向 压力 下 整个 原始 
粗糙 表面 的 真实 接触 面积 。 


1.2 小 尺度 部 分 


在 大 尺度 模型 中 ,每 个 接触 面积 单元 被 认为 是 
完全 接触 的 光滑 平面 , 当 考 虑 短波 长 粗糙 的 影响 时 ， 
上 一 节 得 到 面积 单元 将 不 再 可 以 视 为 完全 接触 。 对 
于 某 一 面积 单元 ,其 真实 接触 面积 是 由 施加 于 其 上 
面 的 法 向 载 丛 与 其 小 尺度 几何 形 貌 决定 的 。 而 其 中 
法 向 载荷 已 经 通过 上 一 步 得 到 。 由 于 粗糙 表面 被 假 
设 为 均匀 的 和 各 向 同性 的 ,因此 可 以 进一步 假设 每 
个 面积 单元 上 的 形 貌 具有 相同 的 特征 , 且 可 以 用 高 
频 PSD 函数 唯一 的 表征 。 借 助 PSD 函数 ,本 研究 将 
每 个 接触 面积 单元 上 的 微 形 貌 与 分 形 粗 燃 表 面 模型 
建立 联系 ,并 进一步 得 到 它 的 力学 响应 , 即 对 应 法 向 
载 集 下 的 真实 接触 面积 。 

正如 式 (1) 所 描述 的 ,频率 在 g, ~ ,的 范围 内 的 
PSD 函数 呈 对 数 线 性 减 小 的 趋势 ,对 应 尺度 的 粗粮 
表面 具有 分 形 的 特征 。 将 PSD 函数 在 9, ~ 9, 的 范围 
内 在 对 数 轴 下 均匀 离散 成 m 部 分 ,每 个 部 分 的 频率 
是 前 一 个 的 万 售 。 对 于 第 i 个 频率 段 的 下 限 9 ,有 

0 = kg, 0<ism-l (2) 
而 大 可 以 由 下 式 得 到 


t= (3) 
gq 


假设 每 一 个 面积 单元 (大 尺度 部 分 给 定 ) 的 形 
貌 可 以 描述 为 m 层 不 同 波长 的 粗糙 峰 的 全 加 。 鉴 
于 此 ,采用 m 阶 二 维 Cantor 集 建 模 表征 在 每 个 面积 
单元 上 其 的 几何 形 貌 。 

如 图 2 所 示 , 上 述 频率 段 所 对 应 的 粗糙 峰 被 转 
化 成 了 二 维 Cantor 集 对 应 的 某 一 层 和 矩形 上 估 。 和 矩形 齿 
块 体 的 边 长 上 是 对 应 的 粗糙 转化 成 的 方 波 的 半 波 
长 ,高 度 h 是 频率 在 g; ~ 4, 范 围 内 的 方 波 按 一 定 方 
式 徐 加 的 等 效 幅 值 之 和 的 2 倍 ,分 别 可 以 表示 为 


1 _ 2 
L 9， 和 
qs 
h, =2 |0.|C(q)gdq (4) 
4 


其 中 ,0. 为 等 效 转 化 系数 ,是 一 
的 量 。 


个 与 表面 参数 有 关 
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间 铎 ， 等 : 粗糙 对 


| 
图 2 二 维 Cantor 集 


Fig.2 Two-dimensional Cantor set 

为 了 确定 9. 的 值 ,首先 要 确定 给 定 频率 段 对 应 
的 方 波 组 合成 Cantor 集 的 麦 加 形式 ,而 这 些 方 波 的 
幅 值 与 PSD 函数 之 间 通 过 一 定 方式 进行 关联 ,从 而 
最 终 决定 了 万 的 值 。 假 设 频率 为 ;的 方 波幅 值 为 
B ,为 了 将 幅 值 为 B;,, 的 方 波 以 保持 9 的 方 波 的 
峰 = 峰 值 不 变 的 条 件 炙 加 在 其 上 , 则 需 额外 状 加 一 
代 来 与 4 相同 ,相位 差 为 且 幅 值 为 B,, 的 广 
波 〇 假设 下 一 层 的 方 波 总 是 关 加 在 上 一 层 的 波峰 范 
甬 内 , 波 谷 处 的 释 加 被 忽略 。 按 这 样 的 方式 逐 层 秋 
加 和 希 同 频率 的 方 波形 成 Cantor 集 , 则 hh. 与 实际 所 舍 


的 直 波 幅 值 5, 可 表示 为 
所 h, = 2B,, 
CN _ h; — hi _ Ah, 


0 = Bb, 一 了 1 三 (5) 
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其 前,A 太 是 相 邻 两 层 Cantor 集 块 体 之 间 的 间距 。 
-三 为 了 让 Cantor 集结 构 与 原始 的 小 尺度 分 形 结构 
应 重大 程度 的 相同 ,考虑 它们 的 PSD 之 间 的 关联 。 
由 学 Cantor 集 模型 中 仅 包含 了 形 如 9, = 4, = 2/24 
的 志波 ,其 PSD 函数 不 再 是 各 向 同性 的 ,因此 令 其 
在 品 gq, = 4, 方向 上 与 原来 的 分 形 PSD 函数 相等 。 忽 
略 相位 的 影响 ,分 解 方 波 5(q) 为 一 系列 的 正弦 波 并 
在 g 方 向 的 坐标 轴 上 秋 加 ,得 到 在 频率 为 g 处 的 正 
弦 波 幅 值 为 


a(g) = Jb(g) (6) 
其 中 ， 
n= 二 (71 (7) 


则 Cantor 集 在 g 处 对 应 的 PSD 函数 值 为 
(| 
(£1) 


Cl(q) = 


_ B87 /h(g) -hl(g +Ag)Y’ 
= A ) (8) 


Ce 面 接触 的 层级 模型 


当 Ad 足够 小 时 令 C,,(g) =C(09) ,将 式 (4) ~ 
(6) 代 入 到 式 (8) 中 ,最 终 得 到 0, 的 值 , 即 
下 
0. = 16H”’ (9) 
由 此 ,Cantor 集 的 几何 形 貌 完全 由 表面 的 小 尺 
度 分 形 特征 确定 下 来 。 
当 法 向 载荷 施加 在 上 述 得 到 的 Cantor 集 模型 上 
时 ,假设 位 于 上 层 的 块 体 先 发 生 接触 变形 ,并且 直到 
位 于 某 一 层 之 上 的 所 有 层 被 完全 压 平 ,这 一 层 才 开 
始 接触 并 发 生 法 向 的 变形 。 随 着 法 向 载荷 的 增加 ， 
Cantor 集 被 逐 层 压 平 ,接触 面积 不 断 增 加 ,直到 最 下 
面 的 一 层 进入 接触 ,这 时 这 个 接触 面积 单元 被 认为 
是 完全 接触 的 。 块 体 在 变形 过 程 中 的 体积 守恒 被 忽 
略 了 。 每 个 特定 的 法 向 载荷 值 对 应 着 Cantor 集 某 一 
被 压 平 的 层 ,也 即 对 应 着 特定 的 接触 面积 。 
第 i 层 Cantor 集 对 应 的 面积 可 以 表示 为 
A = 人 +1 
, 2k 
其 中 ,4 为 最 下 面 那 一 层 的 面积 ,与 对 应 接触 单元 的 
面积 相等 ,为 


2 
) 4， 0<i<sm-1 (10) 


2 
4 = 所 - (11) 
qd 


为 了 得 到 每 层 的 临界 承载 压力 ,一 个 面积 等 效 
的 圆柱 刚性 平 压 头 法 向 压 入 弹性 半 无 限 基 体 上 的 解 
析 解 被 使 用 。 对 于 半径 为 +; 的 圆柱 压 涉 压 入 基 
体 深 度 为 A 访 的 情况 ,平均 接触 压力 Am 表示 为 
Ap; = 2E'Ah,/( mr;), 


(12) 
其 中 := EA/(1 -vw ) 为 等 效 杨 氏 模 量 ;E 为 材料 的 
弹性 模 量 ;z 为 材料 的 泊 松 比 。 
考虑 到 实际 粗 烙 表面 PSD 分 析 得 到 的 频率 是 
连续 的 ,离散 数 m 应 趋 近 于 无 穷 大 。 因 此 ,每 层 的 
接触 面积 4 和 对 应 的 临界 压力 dp 可 表示 为 
gq 
2F’ 
dp = a (13) 
则 压 平 Cantor 集 模型 中 具有 接触 面积 4 对 应 的 
层 的 法 向 临界 总 压力 为 


， 厅 6 
Fe = [Adp = 4B /qr Vm) 9 -qr™ 


(14) 
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结合 在 第 1.1 节 中 得 到 的 每 个 接触 单元 上 法 向 
载 答 ,每 个 接触 单元 的 真实 接触 面积 均 可 以 相应 地 
计算 得 到 。 将 每 个 接触 单元 的 真实 接触 面积 相 加 ， 
即 可 得 到 原始 粗糙 表面 的 总 体 真实 接触 面积 。 


2 有 限 元 模拟 验证 


为 了 验证 上 述 模型 的 准确 性 ,本 研究 利用 有 限 
元 模拟 计算 了 具有 不 同 参数 特性 的 粗糙 表面 的 法 向 
接触 。 利 用 ABAQUS ,构建 了 如 图 3 所 示 的 表面 粗 
糙 的 弹性 基体 模型 。 基 体 的 横向 尺寸 和 表面 的 网 格 
数量 与 粗糙 表面 的 尺寸 与 采样 点 数 一 致 。 粗 糙 表 面 
的 数据 用 修改 节点 坐标 的 方式 加 在 基体 的 上 表面 ， 
并 在 厚度 方向 线性 插值 。 靠 近 表面 的 区 域 由 非 协调 
六 面体 单元 (C3D81) 组 成 。 基 体 侧面 边界 的 节点 沿 
ee ne 
被 施加 在 基体 上 方 一 定 距 离 的 刚性 平面 上 。 接 触 界 
面 酌 切 向 行为 设 为 无 摩擦 条 件 。 弹 性 模 量 和 泊 松 比 
分 别 为 E = 100 GPa,v = 0.3。 


202307 加 


chinaXiv 


图 3 粗糙 表面 的 有 限 元 模型 
Fig.3 FEM model of the rough surface 
在 每 组 模拟 中 ,计算 2 个 粗糙 表面 在 相同 条 件 
下 的 接触 响应 ,其 中 第 1 个 表面 是 PSD 函数 形 如 式 
(1) 的 原始 表面 ,第 2 个 表面 是 用 1.1 节 的 方法 处 理 
得 到 的 截断 表面 。 其 中 ,原始 表面 的 规模 是 256 x 
256 的 采样 点 。 存 有限 元 中 ,为 了 平滑 地 过 渡 表 面 


高 度 ,表面 的 六 面体 线性 插值 网 格 为 10 层 , 整 体 地 
网 格 总 数 接近 10* ,而 经 过 截断 处 理 后 的 网 格 表面 为 
32 x32 或 更 小 (由 截断 频率 而 定 ) ,截断 表面 的 网 格 
总 数 不 足 10 ,计算 规模 减 小 了 10" 倍 以 上 。 将 第 2 
个 表面 上 每 个 接触 节点 的 法 向 力 按 步 提 取 , 然 后 代 
入 到 式 (14) 中 计算 对 应 面积 单元 的 真实 接触 面积 
并 求 和 。 所 得 到 的 力 - 面积 响应 曲线 与 原始 表面 直 
接 用 有 限 元 计算 得 到 的 曲线 进行 比较 ,部 分 结果 如 
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图 4 所 示 。 
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(©) 天 0.8.o=0.5 Lm,g—=5 X 10° rad/m,q.=5.03 X 10° rad/m,qg=4.02 X 10° rad/m。 
妈 4 有 限 元 计算 的 原始 表面 和 截断 表面 
模型 修正 后 的 力 -面积 曲线 
Fig.4 The Force-Area curves of the initial surface 
and the cut-off surface calculated by FEM 

结果 表明 ,对 于 不 同 参数 的 粗糙 表面 ,层级 模型 
预测 的 力 - 面积 曲线 与 原始 粗糙 表面 的 响应 结果 相 
较 于 其 他 理论 模型 更 加 接近 。 这 说 明 本 研究 提出 的 
模型 在 预测 多 尺度 粗 烟 表面 的 法 向 接触 响应 中 是 可 
行 的 。 在 实际 应 用 中 ,针对 多 斥 度 粗糙 的 各 向 同性 
表面 ,只 需 提取 一 定 精度 的 信息 得 到 较 大 尺度 的 轮 
廓 ,以 及 小 尺度 的 分 形 特征 ,利用 本 研究 模型 即 可 在 
较 小 的 计算 成 本 下 获得 其 法 向 接触 的 响应 。 该 模型 
可 以 通过 改变 材料 属性 扩展 到 弹 塑性 情况 ,这 一 部 
分 将 在 之 后 的 工作 中 进行 分 析 。 
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